Ising phases of Heisenberg ladders in a magnetic field 
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We show that Dzyaloshinskii-Moriya (DM) interactions can substantially modify the phase dia- 
gram of spin-f/2 Heisenberg ladders in a magnetic field provided they compete with exchange. For 
non frustrated ladders, they induce a local magnetization along the DM vector that turns the gap- 
less intermediate phase into an Ising phase with broken translational symmetry, while for frustrated 
ladders, they extend the Ising order of the half-integer plateau to the surrounding gapless phases of 
the purely Heisenberg case. Implications for experimental ladder and dimer systems are discussed. 

PACS numbers: 75.10.-b, 75.10.Jm 



Spin ladders [l[ have played an important role in the 
field of strongly correlated systems over the past fifteen 
years. They are among the best studied spin-gap sys- 
tems, thanks to the numerous experimental realizations 
in transition metal oxides 0] and organometallic com- 
pounds their intermediate phase in a magnetic field 
is one of the simplest realizations of a Luttingcr liquidd, 
and frustrating them with diagonal rungs opens a mag- 
netization plateau at half-saturationQja simple example 
of non- integer magnetization plateaux[6(. 

All this is true for purely SU(2) Heisenberg interac- 
tions, but in actual systems, anisotropic terms such as 
the Dzyaloshinskii-Moriya (DM) interaction are often 
present. In ordered antiferromagnets, they open a small 
gap in the spin-wave spectrum and may induce a small 
canting, but otherwise they can be neglected for most 
purposes. One may thus be tempted to conjecture that, 
for gapped antiferromagnets, their effect will be negligi- 
ble except in the gapless phase induced by a magnetic 
field, where they will open a gap. The intriguing prop- 
erties of the famous plateau system SrCu2(B03)2|7|, in 
particular the recent report that the translational sym- 
metry remains broken above the 1/8 plateau [H, suggest 
this might not be so simple however. Indeed, significant 
DM interactions have been identified in this system, and 
the role they may play between the plateaux remains to 
be investigated. In the same spirit, the nature of the 
phase transition of the intermediate phase of the spin 
gap system (^^(CsH^^^CLiQ is still an open issue, 
and the recent identification of DM interactions in that 
system fioj] opens new perspectives that remain to be ex- 
plored. 

In this Letter, we report on an extensive investigation 
of the effect of DM interactions in gapped systems in the 
context of spin ladders, and we show that even very small 
DM interactions can induce Ising phase transitions with 
broken translational symmetry in all spin ladders, frus- 
trated or not, under appropriate conditions. To cover all 



FIG. 1: Sketch of the frustrated ladder of Eqs. 0-©. The 
thick solid lines stand for J, the thin ones for Jy , the dashed 
ones for Jx, and the arrows for D. We have also shown the 3 
generators of the space group: the a xz and o y ~ mirror planes 
and a twofold rotation axis Ci(z). 



aspects of the problem, we concentrate on a fairly general 
model depicted in Fig. [Tj Its Hamiltonian can be split 
into single-rung exchange, DM, and Zeeman energies, 
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and Heisenberg interactions between adjacent rungs, 

'Hint = J\\ • + + Sj,2 ' Sj + i,2) 
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+Jx / aSj-.i ■ Sj+1,2 + Sj,2 ■ Sj'+i.i) (2) 
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In this work, the magnetic field H = (0, 0, H) is applied 
perpendicular to the DM vector D = (0,D,0), which is 
taken to be uniform, the relevant configuration for lad- 
ders that are translationally invariant but have no inver- 
sion center in the middles of the rungs. The full sym- 
metry group of the system can be generated from three 
symmetry operations: Ta yz , Ci(z) and Ta xz , as shown 
in Fig. [T] where the additional time-reversal operator T 
is needed to compensate for the presence of the magnetic 
field. The symmetry operations act both in spin and real 
space. 
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The physics of this model is controlled by: i) The level 
of frustration of the Heisenbcrg exchange. Indeed, when 
D = 0, a magnetization plateau at half the saturation has 
been shown to open in strongly frustrated ladders (1/3 < 
Jx/J\\ < 3), corresponding to a gapped phase with a 
broken translational symmetry^. By contrast, non frus- 
trated ladders (Jx = or Jm = 0) and weakly frustrated 
ladders (0 < Jx/J\\ < 1/3 or < J\\/Jx < 1/3) remain in 
a gapless Luttinger liquid phase from zero magnetization 
to full polarization, ii) The possible competition between 
exchange and DM interactions fill] . In the present con- 
text, this can be understood as follows. Due to the DM 
interaction, an isolated rung in a field develops opposite 
moments along D x H, on the two sites of the dimerflil]. 
If such dimers are coupled in a ladder geometry, an an- 
tiferromagnetic inter-rung coupling J|| will compete with 
the DM interaction because Jii tends to align the mo- 
ments antiparallel along a leg. By contrast, an antiferro- 
magnetic inter-rung coupling Jx will not compete since 
both Jx and the DM interactions induce antiparallel mo- 
ments along the diagonals. 

The non-competing case (Jx/J|| > 1) has been stud- 
ied previously |l3j. The DM interaction opens a gap in 
the gapless incommensurate phase, as expected, but only 
induces a small shift of the Ising plateau transition in 
strongly frustrated ladders. In the following, we will con- 
centrate on the competing case (Jx/J|| < !)■ 

The simplest way to map out the phase diagram is to 
use a Hartree variational function of the form 



^ = Jj(cos||T) 



sin ^ |S> 



(3) 



with \T) = | hU) and \S) = (| I1T2) - | Tii 2 »/v^ on 
the j-th rung [14j . The trial state is characterized by the 
symmetric and antisymmetric magnetizations of a rung, 
m j± = ± Sj,2). The nonzero components are 



irrij _ = sin t 



e lip j/V2, m z jt+ =cos 2 (6» i /2). (4) 



The variational phase diagram without DM interac- 
tion is shown in Fig. [Ha). It is symmetric with re- 
spect to the exchange of Jx and J|| and contains all the 
phases found previously [13| : two fully symmetric inte- 
ger plateau states with m| + = and 1; a half-integer 
plateau phase with a broken Z2 translational symmetry, 
rrij + — 1/2 cx (— and a gapless phase with a sponta- 
neously broken 0(2) symmetry (the latter being an arti- 
fact of the mean- field approximation). 

The variational phase diagram in the presence of the 
DM interaction is shown in Fig.[5{b). It is clearly related 
to the 0(2)-symmetric case: some of the transition lines 
lie close to their D = counterparts. However, there are 
some notable differences. In the noncompeting regime, 
Jx/J\\ > lj the phase transitions between the integer 
plateaux and the intermediate gapless phase are re plac ed 
with crossovers, in agreement with previous work 




FIG. 2: (Color online) Variational phase diagram for (a) D = 
and (b) D/(Jx + J\\) = 0.1. On the left and bottom axes, 
Jx and H — Ho, where Ho = J + (Jx + J||)/2, are given in 
the units of Jx + J11. In (a), the spins forming singlets and 
triplets in the gapped phase are represented by open and solid 
circles. In the gapless 0(2) phase the arrows shows a typical 
mean-field ground state. In (b), the 0(2) symmetry is broken 
by finite D: a unique state is selected in the non frustrated 
case, an a two-fold degenerate state in the frustrated case. 
The Ising phase of the plateau in (a) develops continuously 
by turning on D. Along the red dashed line the site-factorized 
wave function is an exact wave function in the large J limit. 
R and R' denote the tricritical points. 



Far more dramatic changes can be seen in the compet- 
ing regime, Jx/J\\ < 1- At strong frustration, 1/3 < 
Jx/J\\ < 1, the Ising-ordered phase has substantially 
expanded its boundaries by absorbing the incommen- 
surate phases. A distinct Ising order appears for zero 
or weak frustration, < Jx/ J\\ < 1/3. In the upper 
Ising phase, the order parameter is the staggered com- 
ponent of the magnetizations m| + and nij _ . It breaks 
Ta xz . In the lower Ising phase, the order parameter is 
the staggered component of rrij _ . It breaks Ta yz . Since 
T <J yz Ci(z)T o xz is the translation operator, both phases 
break the translational invariancc as well. In any of the 
two phases, only one generator of the space group is bro- 
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FIG. 3: (Color online) (a) Uniform magnetizations of neigh- 
boring rungs (mo+ and mf + ) and average (m) in the upper 
Ising phase and (b) staggered magnetizations m§ + and mi + . 
Solid lines: variational approach; symbols and dotted lines: 
DMRG for 58 sites; shaded area: plateau phase when D = 0. 
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FIG. 5: (Color online) Schematic representation of the phase 
transitions. Away from the A = 1 Heisenberg line, the 7i e s 
loses the axial symmetry around the total field: the spins 
develop staggered order along the long dashed blue arrows if 
A < 1 and along the short dashed red arrows if A > 1. 




FIG. 4: (Color online) Magnetizations in the lower Ising 
phase. Solid lines: variational approach; symbols and dot- 
ted lines: DMRG for 80 sites. 



system to maintain a non-zero staggered magnetization 
along the D vector. We find a qualitative agreement 
between the variational and numerical results. Quanti- 
tatively, the mean-field approach overestimates the order 
parameter by as much as a factor of 2, presumably the 
result of neglecting quantum fluctuations. 

To understand the physics of the ordered states and 
phase transitions between them we turn to the strong- 
coupling limit where the rung exchange J of Eq. ([lj) dom- 
inates. The low-energy physics of the model then reduces 
to that of a spin-1/2 XXZ chain 



K 



efi 



(5) 



ken, thus the transitions from the uniform phase into 
these phases are continuous. At the variational level, 
the transition between the two Ising phases is first order, 
but numerical results and field-theory arguments suggest 
that the transition is made continuous by quantum crit- 
ical fluctuations (see below). 

To test the reliability of the variational approach, we 
have obtained the ground state and magnetization dis- 
tributions for the ladder model (p][2|) using the Density 
Matrix Renormalization Group (DMRG) for two repre- 
sentative cases: the unfrustrated ladder (Jx = 0) and 
a strongly frustrated one (Jx/J\\ = 2/3), and we have 
obtained strong evidence for the existence of the two 
Ising phases. The results for relevant local magnetiza- 
tions are compared with the variational results in Figs. [3] 
and[U In the upper Ising phase (Jx/J\\ = 2/3, Fig. [3]), 
the agreement is nearly perfect and the presence of a 
broken symmetry far outside the half-integer magnetiza- 
tion plateau is well established: magnetizations m| + of 
two neighboring rungs exhibit a finite difference in al- 
most the entire range of fields between the two integer 
magnetization plateaus. For Jx/ J\\ = (sec Fig. @|, we 
have included a small staggered field to allow the finite 



with j = J|| - Jx and A = ( J x + J|j)/2(J|| - Jx)- The 
longitudinal field is related to the field in the original 
problem, h z = H — J — (Jx + J\\ )/2, while the transverse 
field comes from the DM interactions, h x = D/V2, as 
discussed previously [Hj]. The spin-1/2 operators a arc 
related to the original spin operators in the following way: 



P 



y/2 



where a = x, y and P = (S^ - 5+ 2 )(Sr 1 - Sj 2 )/2 is the 
projector onto the subspace \T)}. Nonzero averages 
(a?) and (ct^) signal the appearance of antisymmetric 



magnetizations _ and 



while (a z ) translates into 



the symmetric magnetization tti* + — 1/2, cf. Eq. ((4]). 

The XXZ model in a tilted uniform field ([5|) has been 
studied previously [ljl, 17 1. The magnetic order is lost 
across the line of Ising transitions, h z w (A + into 
the fully magnetized plateaus. The two Ising phases are 
separated by the Heisenberg line A = 1 along which the 
system has an axial symmetry. For h = |h| < 2j, the 
system is gapless and can be viewed as a Luttinger liq- 
uid of magnons, whose dimcnsionlcss compressibility K 
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[J] varies between 1 (a dilute magnon gas near the fully 
magnetized state) and 1/2 (a dense magnon fluid at low 
magnetization) [1 81 ] - 

To discuss the physics away from the Heisenberg point, 
it is convenient to make a global rotation in the xz plane, 
defining a new quantization axis z 1 aligned with the total 
field h. When A / 1, this rotation generates new terms 
in the Hamiltonian. The nature of the phase transitions 
becomes particularly transparent in the vicinity of the 
tricritical points R and R' located at A = 1, h = ±2j, 
where the system can be described as a dilute magnon 
gas. Then the main effect is the lowering of the axial 
symmetry by a weak anisotropy induced in the x'y plane, 



i 



A 



+i 



-i). 



. (6) 

where A ~ (A — l)h x 2 /8j. The magnon Luttingcr liquid 
becomes gapped by developing a stag gere d magnetization 
along y if A < 1 (breaking Ta xy [19() or along x' if 
A > 1 (breaking Ta xz [HI), see Fig. El The anisotropy 
([6]) represents a pairing field for the magnons. From the 
scaling dimension of the pairing operator we infer that 
the energy gap in the magnon spectrum should scale as 
Eg oc with a nonuniversal exponent v = K/(2K— 1). 
The global rotation also generates a 3-body term [lq . 



(A - l)h x 



which has only a minor effect on the physics near the crit- 
ical line separating the Ising phases. While low-energy 
magnons live near lattice momentum tt, this term cre- 
ates 1- and 3-magnon states with a total momentum 0, 
which are thus high-energy excitations. Therefore the 

3- body term does not affect the low-energy physics to 
the first order. In the second order it generates 2- and 

4- magnon terms that simply renormalize the magnon ve- 
locity, compressibility K, and pairing field A but do not 
change the physics qualitatively. The existence of two 
distinct ordered phases is confirmed by an exact solu- 
tion on a special line in the phase diagram where the 
mean- field approach yields the correct result [l6j] . The 
Heisenberg line A = 1 is guaranteed to be critical by the 
0(2) symmetry. 

To summarize, we have shown that even a weak DM 
interaction can substantially alter the phase diagram of 
ladders, provided that it competes with exchange. For 
strongly frustrated ladders, the phase with a broken 
translational symmetry expands to include the incom- 
mensurate regions surrounding the fractional plateau, 
while for unfrustrated and weakly frustrated ladders, the 
gapless phase undergoes an Ising transition that also 
breaks translation, but with a different symmetry. Be- 
yond ladders, we expect to find similar effects in higher 
dimensional models of coupled spin- 1/2 dimers, whenever 



DM and exchange compete. In fact, the broken transla- 
tional symmetry above the 1/8 plateau of SrCu2(B03)2 
might be an example of the frustrated case, and the or- 
dered intermediate phase of Cu2(C5Hi2N2)2Ci4 of the 
unfrustrated one. However, taking into account the spe- 
cific (and complex) geometries of these systems is neces- 
sary to go beyond these qualitative observations. 
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